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ABSTRACT 

Aims. We aim to study the properties of the dense molecular gas towards the inner few 100 pc of four nearby starburst galaxies dom- 
inated both by photo dissociation regions (M82) and large-scale shocks (NGC253, IC342 and Maffei2), and to relate the chemical 
and physical properties of the molecular clouds with the evolutionary stage of the nuclear starbursts. 

Methods. We have carried out multi-transitional observations and analyses of three dense gas molecular tracers, CS, HC 3 N 
(cyanoacetylene) and CH3CCH (methyl acetylene), using Boltzmann diagrams in order to determine the rotational temperatures and 
column densities of the dense gas, and using a Large Velocity Gradients model to calculate the H 2 density structure in the molecular 
clouds. 

Results. The CS and HC3N data indicate the presence of density gradients in the molecular clouds. These two molecules show similar 
excitation conditions, suggesting that they arise from the same gas components. In M82, CH3CCH has the highest fractional abun- 
dance determined in a extragalactic source (1.1 x 10~ 8 ). 

Conclusions. The density and the chemical gradients we have found in all galaxies can be explained in the framework of the starburst 
evolution, which affects the chemistry and the structure of molecular clouds around the galactic nuclei. The young shock-dominated 
starburst galaxies, like presumably Maffei 2, show a cloud structure with a rather uniform density and chemical composition which 
suggests low star formation activity. Molecular clouds in galaxies with starburst in an intermediate stage of evolution, such as NGC 253 
and IC 342, show clouds with a large density contrast (two orders of magnitude) between the denser regions (cores) and the less dense 
regions (halos) of the molecular clouds and relatively constant chemical abundance. Finally, the galaxy with the most evolved starburst, 
M 82, has clouds with a rather uniform density structure, large envelopes of atomic/molecular gas subjected to UV photodissociating 
radiation from young star clusters, and very different chemical abundances of HC3N and CH3CCH. 

Key words, galaxies: starburst - ISM: molecules - galaxies: individual: NGC 253, M 82, IC 342, Maffei 2- galaxies: nuclei - Radio 
lines: galaxies 

1 . Introduction the molecular clouds. Such studies also provided a powerful tool 

to study not only the type of dominant activity but the evolution 
Molecular clouds within the central region of our galaxy pro- of the dense ISM in the obscured regions within galactic nu _ 

£ ■ vide good templates for a better understanding of the molec- dei For instanc6i the H CO + -to-HCN line intensity ratios have 
CG ; ular emission from the nuclei of external galaxies. The study been used for discr i m i nating between AGNs and starburst sig- 
of well known molecular clouds in the Milky Way, like those natur es in the galaxy ce nters dKohno et al]l200lr. IrrrTanishi et al. 
affected by shock waves due to cloud-cloud collisions, mass 2 004: lKrins et al. 2008). Other intensity ratios, like HNC/HCN, 
loss from massive stars, gas accretion and explosive events like ^ also used for d i fferent i ating between XDRs and PDR s con- 
su£ernovae (e.g. | Martm-Pintado et al. || 1997t | Htittemeister et al. | trihuiion in different stages of evolution (Loenenet al. 2008; 
L1998J), those pervaded by UV dissociating radiation from the Baan et al 2008) As another example which uses column den . 
massive stars that create large photo dissociation regions (PDRs) sities derived from multi _ line ana i yses , the HNCO/CS abun- 
(e.g.| Martfnetal. || 2008a|) , and also those which study regions dance ratio was proposed for differentiating betwee n the star- 
affected by X-rays (XDRs) (e.g. | Martfn-Pintado et al. | |200fi bursts maMy dominated by shocks or UV fields. (iMartm et all 
lAmo-Baladron et al.H2009t) . allows to use the chemical compo- 2009a) 
sition to infer the physical processes dominating the heating of 

" ' The nearby (~3Mpc) galaxies NGC 253 and M82, are 

a of the brightest infrared extragalactic sources, and among 
the most outstanding in terms of detec t ion of extragalactic 
molecules (Mauersberg er & Henkelll993t IMartm et alJ l2006b. 



* This work is based on observations with the IRAM 30-m telescope, two of the brightest infrared extragalactic sources, and among 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and the most outstanding in terms of detec t ion of extragalactic 
IGN (Spain). 



2 



Aladro et al. (2010): CS, HC 3 N and CH 3 CCH towards the nuclei of starburst galaxies 



Aladro et al., 2010b, in prep.). Simple molecules such as CO or 
CS are observed towards both sources with similar column den- 
sities and abundances. However, nearly all the species detected 
so far, in particular the more complex ones (NH3, HNCO or 
CH3OH), as well as SiO, show systematically lower abundances 
in M 82 than in NGC 253 and other starburst galaxies such a s 
IC342 and Maffei2 (iTakano et al.l 12003b IMartfn et al.ll2009al) . 
The observed differences among galaxies are not only due to ex- 
citation, but reflect different chemical compositions. The large 
overall abundance of HCO in M 82 has been interpreted as ev- 
idence for a chemistry mainly dominated b y photon-dominated 
region s (PDRs) within its central ~650pc (IGarcfa-Burillo et ail 
2002), which does not mean that other proce sses like shocks o r 
X-ray heating do not take place in its nucleus (Baanetal . 2008). 
The dissociating radiation explains the low a bundance in M 82 of 
the complex and fragile spe cies cited above ( Mauersb erger et all 
120031; lMartfn etal]r2006al) . Other molecules like CS are not 
under-abundant in M82 relative to NGC 253, which can also 
be explained in terms of PDRs chemistry, since a high abun- 
dance of this molecule is predicted for slightly shielded regions 
dDrdla et al.lll989t ISternberg & Dalgarnoiri995l) . This evidence, 
as well as the detection of other PDR enhanced species such as 
HOC + and CO + has led to think of M 82 as the pro t otype of a 
giant extragalactic PDR (iNguven-O-Rieu et al.|[l 989; Mao et al. 
|200(* IGarcfa-Burillo et al]l2002t IMartfn et alJl2009ai) . 

On the other hand, the central regions of NGC 253 and 
Maffei2 appear to resemble quite closely the giant molecu- 
lar cloud complexes within our Galactic central regions where 
the low velocity shocks and massive sta r formation dominat e 
the chemistry of the molecular material (Martin et al. 2006b). 
Through the HNCO/CS ratio it has been shown that, the nuclear 
region of NGC253, though less dominated by UV radiation than 
in M82, is significantly more UV pervaded than in IC342 or 
Maffei 2. This has been confirmed by the recent de tection of the 
PDR tr acers HCO, HOC + , and CO + in NGC 253 ( IMartfn et al.l 
2009b). The observed chemical differences are interpreted as the 
evolution of the starburst phenomenon, where M 82 would be in 
a later stage (probably in a post-starburst phase where almost all 
the gas has already been converted into stars), while NGC 253 
or Maffei 2 centers would host younger starbursts (where still 
there are a large amount of reservoir gas which is being used to 
form the stars). This point can be used to link the stage of evo- 
lution of the starburst galaxy centers with the dominant physical 
processes (PDRs, shocks, XDRs), in such a way that a possible 
scenario, as seen by the single dish telescopes, could be that the 
more evolved (or post-starbursts) galaxies, like M82, are mainly 
dominated by UV fields creating large PDRs, while the younger 
starburst are basically dominated by shocks between the molecu- 
lar clouds. It does not necessary implies that these are the unique 
processes taking place, but the dominant ones. 

Interferometric maps of several molecules towards the center 
of IC 342 show that both, PDRs and shock-do minated regions, 
can b e resolved within its inner few lOOpc (Mei er & Turner! 
120051) . while observations with single dish telescopes do not 
have the resolution to distinguish between these two types 
of activities in the nucleus. Thus, the overall scenario for 
IC342, based on beam averaged observations, is a mainly 
shock-dominated nu cleus, as shown by the HNCO/CS ratio 
(Mar tin et al.ll2009ab . This galaxy seems to be in an early star- 
burst sta ge, as shown by some shocks tracers like SiO or 
CH3OH (lUsero et al.ll2006l) . On the other hand, Maffei 2 seems 
the galaxy in our sample where C-shocks produced by cloud- 
cloud collisions play the most important role, being able to ex- 



plain its high gas kinetic temperatures and chemical abundances 
(iMauersberger et al .1120031) . 

In this paper we present observations of CS, HC3N and 
CH3CCH towards the nuclei of four starburst galaxies. We aim 
to relate the chemical and physical properties derived from the 
analyses of these three dense gas tracers (densities, tempera- 
tures and fractional abundances) with the stage of evolution of 
each starburst. The outline of the paper is as follows. In Sect. [2] 
we present our observations and other complementary line data 
taken from the literature. In Sect. [3] we explain the data analy- 
sis. The data are first analyzed under the local thermodynamic 
equilibrium (LTE) approximation in Sect. 13.11 where we obtain 
the column densities and rotational temperatures, which give us 
a first hint on the cloud structure. Sect. I3.2| presents the results 
from a LVG modelling under non-LTE assumption. From this 
model, we estimate the H2 densities for different cloud compo- 
nents. In Sect. |U we discuss the results of both approaches, tak- 
ing into account their validity and limitations, and relating those 
results with the evolutionary stage of the molecular clouds in 
each galaxy. Finally, in Sect.[5]we present our conclusions. 

2. Observations and data reduction 

We have detected a total of 37 transitions of CS, CH3CCH and 
HC3N towards the nuclear regions of the starburst galaxies M 82, 
NGC 253, IC 342 and Maffei 2, plus 5 additional non-detections. 
26 out of the 37 transitions are newly detected. Table Q] shows 
the main characteristics of these galaxies and Table [2] lists the 
observed lines. We have used the IRAM 30-m telescope (Pico 
Veleta, Spain), except for the CH3CCH (13 - 12) line, which has 
been observed with the JCMT telescope (Mauna Kea, Hawaii, 
USAfl The observations were carried out in several periods be- 
tween 2005 and 2009. At the IRAM 30-m telescope, we used 
the, now decommissioned, ABCD SIS receivers tuned simulta- 
neously in single sideband mode in the 1, 2 and 3 mm bands. 
Observations were carried out in wobbler switching mode with 
a symmetrical throw of 220" in azimuth and a switching fre- 
quency of 0.5 Hz. As back-end, we used the 256 x 4 MHz filter 
banks. The beam sizes for each frequency are given in Table [2] 
The pointing was checked every one or two hours towards sev- 
eral standard pointing calibrators, with an accuracy of ~ 3". 
The spectra were calibrated with a standard dual load system. 
The image sideband rejection at the observed frequencies ranged 
from 10 dB to 32 dB. The JCMT observations were carried 
out during the summer period in 2005. Observations were per- 
formed in beam switched mode with a frequency of 1 Hz and a 
beam throw of 2' in azimuth. We used the A3 receiver in dou- 
ble side band to observe the CH3CCH 7=13-12 transition at 
222.1 GHz. As spectrometer we used the now decommissioned 
Digital Autocorrelation Spectrometer (DAS) in wide band mode 
which provided a bandwidth of 1800 MHz at a 1.5 MHz spectral 
resolution. The telescope beam size was 22". The pointing was 
also checked every hour, with an accuracy of ~3". 

The observed positions of the galaxies are shown in Table Q] 
M 82 observations were pointed towards the North-East molec- 
ular clump, where the photodissociating radiation is claimed 
to be stronger, as seen in the HCO interferometric maps of 
IGarcfa-Burillo et all (|2002|) (offsets (+13", +7.5") with respect 
to its coordinates in Table [TJ. In some cases, emission from the 

1 The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organization for 
Scientific Research, and the National Research Council of Canada. 
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Table 1. Main properties of the galaxies in our study. 



Galaxy 


RA 


DEC 


Type 


D a 


Vlsr" 




N(H 2 ) d 


Lir e 


SFR' 


Mm c 


M 8 




(J2000) 


(J2000) 




(Mpc) 


(kms -1 ) 


(") 


(xlO 22 crrT 2 ) 


(xlO 9 Lq) 


(M yr-') 


(xlO 9 M ) 


(xlO 8 M Q ) 


M82 


09:55:51.9 


+69:40:47.1 


10 


3.6 


300 


12.0 


7.9 


29.74 


~9 


1.30 


5.2 


NGC 253 


00:47:33.4 


-25:17:23.0 


SAB(s)c 


3.9 


250 


20.0 


6.2 


15.10 


3.6 


2.57 


28.0 


IC 342 


03:46:48.5 


+68:05:46.0 


SAB(rs)cd 


3.3 


31 


10.6 


5.8 


2.26 


2.5 


18.2 


1.4 


Maffei 2 


02:41:55.1 


+59:36:15.0 


SAB(rs)bc 


3.3 


-17 


11.0 


4.4 


2.7 


0.5 




8.0 



a Values from lFreedman et all 119941) for M 82jKarachents ev et alj 120031) for NGC 253 and IC 342 and lFineerhuTe t al. 1 2007) for Maffei 2. 

b Velocities respect to the Local Standard of Rest. In M 82 the value refers to the North-East lobe. 

c Sizes of the emitting region, taken from: lMartfn et all <2006ah for M 82; This work for NGC 253: lBavet et alj 120061) for IC 342: lMauersberger et alj 119891) for 
Maffei 2. 

d Values taken from Mauersberger et al. ( 2003), and further corrected by the source size in column 7. 

c Values taken from de Vaucouleurs et alj 119911) . 

f Values taken from lstrickland et alj (2004) for NGC 253 and M 82: IWalker et"ai] ff992l) for IC342. 

g Estimation of the global molecular gas mass contained in the source. We have assumed that the gas is in a virial equilibrium. For the estimations we used the 
source sizes S of column 7 and an average linewidths calculated from those presented in Table lA. F1 



Table 2. Observed Lines. 



Transition 


Frequency 


reti/O e g 


n 


n ° 

"crit 


Galaxy 




(GHz) 




V ) 


fern" 3 } 




CS(2- 1) 


97.981 


1.25 


25.5 


1.2 xlO 5 


NGC 253, IC 342 


CS (3 - 2) 


146.969 


1.35 


16.8 


2.0 x 10 6 


IC 342, M 82 


CS(5 - 4) 


244.936 


1.79 


10.1 


4.6 x 10 7 


NGC 253, IC 342, Maffei 2, M 82 


HC 3 N (9 - 8) 


81.881 


1.22 


30.2 


1.6 x 10 b 


IC 342, Maffei 2 


HC 3 N(10-9) 


90.979 


1.23 


27.5 


3.6 x 10 6 


IC 342 


HC 3 N(11 - 10) 


100.076 


1.25 


24.9 


4.8 x 10 6 


IC342 


HC 3 N(12- 11) 


109.174 


1.27 


22.2 


7.2 x 10 6 


IC 342, Maffei 2 


HC 3 N(15 - 14) 


136.464 


1.32 


18.2 


2.6 x 10 7 


IC 342, M 82, Maffei 2 


HC 3 N(16- 15) 


145.561 


1.35 


16.9 


7.0 x 10 7 


IC 342, M 82 


HC 3 N(17- 16) 


154.657 


1.38 


16.0 


1.0 x 10 8 


IC342,M82 


HC 3 N(18- 17) 


163.753 


1.41 


15.1 


1.3 x 10 8 


M82 


HC 3 N (23 - 22) 


209.230 


1.60 


12.0 




Maffei 2* 


HC 3 N (24 - 23) 


218.325 


1.66 


11.5 




IC342,M82 


HC 3 N (28 - 27) 


254.699 


1.92 


9.7 




M82, Maffei 2* 


CH 3 CCH(5 -4 ) 


85.457 


1.22 


29.1 


3.7 x 10 4 


IC 342, Maffei 2 


CH 3 CCH(6 - 5 ) 


102.548 


1.28 


24.2 


1.3 x 10 5 


IC 342, Maffei 2 


CH 3 CCH(8 -7 ) 


136.728 


1.32 


18.2 


6.8 x 10 5 


M 82, IC 342, Maffei 2 


CH 3 CCH(9o-8 ) 


153.817 


1.38 


16.1 


1.1 x 10 6 


M82 


CH 3 CCH(10 - 9 ) 


170.906 


1.44 


14.4 


3.1 x 10 6 


M82 


CH 3 CCH(13 - 12n) 


222.167 


1.25/1.68* 


21.7/11.3* 


1.3 x 10 7 


NGC 253, Maffei 2* 


CH 3 CCH(14 - 13 ) 


239.252 


1.80 


10.3 


2.9 x 10 7 


Maffei 2* 


CH 3 CCH(15 - 14 ) 


256.337 


1.87 


9.6 


4.4 x 10 7 


M82 


CH 3 CCH(16 - 15 ) 


273.420 


2.02 


9.1 


9.4 x 10 7 


M82 


c-C 3 H 2 (2,, 2 - l ai ) 


85.339 


1.22 


29.2 




Maffei 2 


c-C 3 H 2 (5 2 , 4 -4,, 3 ) 


218.160 


1.66 


11.5 




M82 


H 2 CO(2o, 2 -lo,i) 


145.603 


1.35 


16.9 


2.3 x 10 6 


M82, IC 342 


H,CO(3o, 3 -2 , 2 ) 


218.222 


1.66 


11.5 


9.4 x 10 6 


M82, IC 342 


H,CO (3 22 - 2 2 i) 


218.476 


1.53 


11.5 


5.2 x 10 6 


M82 



Last column indicates the galaxies where each transition was observed. Name in italics means a new detection. 
* Name with * means an upper limit of the line intensity. 

11 Critical densities corresponding to a = 90 K. Due to a lack of collisional coefficients in our LVG model for CH 3 CCH, the n clit for this 

molecule were calculated using the collisional coefficients of CH 3 CN, which are very similar. 
b The first value corresponds to NGC 253 observed with the JCMT telescope, and the second value corresponds to Maffei 2 observed with 

IRAM 30-m telescope. 



center of the galaxy was picked up by the beam, in particular for 
the lower frequency transitions (bigger beams) where the spectra 
show a bump at the left of the lines. NGC 253 was observed to- 
wards the same position where the 2 mm frequency survey was 



carried out by iMartm et al.1 d2006bl) . Maffei 2 and IC342 were 
observed towards their centers. 

The data, shown in Figs.Q~Jto [4] are given in a Tmb scale, 
obtained from T* A as Ty^-iF^IB^) T* A . The conversion fac- 
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1 p/A 


IC3N (16-15) 











R HC3N (17-16) 









1 




HC3JN (24-23 ) 










JT1 CH3CCH(9-8) 







CH3CCH(10-9) 





Fig.l. Gaussian profiles fits to observed lines in M82. Fi §- 3 - Gaussian profiles fits to observed lines in IC342. 

Temperatures are in r MB (K). Velocities are respect to the Local Temperatures are in T MB (K). Velocities are respect to the Local 

Standard of Rest. Some of the transitions show a bump located at Standard of Rest - 
the left of the line coming from the emission of the galaxy center. 

This feature is marked in the spectra with the word "CENTER". r 




Fig. 2. Gaussian profiles fits to observed lines in NGC253. 
Temperatures are in 72mb (K). Velocities are respect to the Local 
Standard of Rest. 



tor (F e ff/S e ff) is given in Table [2] for each transition. Baselines 
of order or 1 were subtracted in most cases. The rms of the 
residuals after the baseline subtraction are shown in Table I A. fl 
Gaussians profiles were fit to all detected lines (see Figs.[TJto |4|. 
We have used the CLASS0 software package for data reduction 
and Gaussian profile fitting. 

The intensity of the molecular emission for all species was 
also corrected by the beam dilution effect, due to the coupling 
between the source and the telescope beam, as 72b = [(0^ + 
(%)/(%] Tmb, where T% is the source averaged brightness tem- 
perature, # s is the source size, Ob is the beam size in arc sec- 
onds, and Tmb is the measured main beam temperature. Since 
CS, HC3N and CH3CCH require similar excitation conditions 
(see critical densities in Table 0, it is plausible to use the same 
source size for all transitions. In the case of NGC253, we com- 
pared the bright ness temperature of CS(2 - 1) observed with the 
SEST telescope dMartfn et al.l2005l) with the same line observed 
with the IRAM 30-m (this work). Using the relation between 
Tmb and 7b, the comparison of both observations leads to a 
20" equivalent sourc e size, which agrees with the value used by 
Mart in et al.l (l2006bl) . For the other galaxies, there are no tran- 
sitions observed at the same position with different telescopes, 




CLASS http://www.iram.fr/IRAMFR/GILDAS 



Fig. 4. Gaussian profiles fits to observed lines in Maffei2. 
Temperatures are in Tmb (K). Velocities are respect to the Local 
Standard of Rest. 



so we adopted source sizes derived from other studies listed in 
Table [TJ The influence of the source size in the rotational tem- 
peratures and column densities can be seen in Ba vet et ai] (f2009) 
for CS in M 82, NGC 253 and IC 342, among other galaxies. The 
assumption of similar source size for all transitions might have a 
small contribution to the derived temperature structure (< 30%) 
and less than a factor of 2 in the total column densities. On the 
other hand the assumed source size for each galaxy will have an 
almost insignificant effect on the temperatures (< 7%) but might 
have a strong effect in the column densities in the case of the 
source size being much smaller {6 S < 2") than the one assumed. 

Apart of CS, HC3N and CH3CCH, other molecules are de- 
tected in the observed bands. In particular, H2CO(2o ; 2 - lo,i), 
H2CO(3o,3 - 2o,2 and H2CO(32,2 - 22, 1) have been detected in 
M82 and IC 342. c-C 3 H 2 (5 2 ,4 - 4 U ) is seen in M82, and c- 
C3H2 (2 12 - 1 1,0) is detected in Maffei2. These lines are shown 
in Figs.[TJto Hand listed in Table[2] No lines are blended except 
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for HC 3 N(16 - 15), which is contaminated by H 2 CO (2 , 2 - ln,i). 
In this case, we estimated the contribution of H2CO (2n,2 - lo,i) 
from other detections of this molecule (Aladro et al. 2010b, in 
prep.) and subtracted it. Then, we fitted a Gaussian profile to the 
residuals. 

In our analysis we also included some results from the 
literature in the 3, 2 and 1 mm bands. We also included the 
CS(7 - 6) sub-millimeter line when available. Details can be 
foun d in Table lATl For NGC 253, we took the HC 3 N lines 
from Mauersberger et all (Il990h . that cover the full range from 
81.9 GHz to 236.5 GHz. W e also used several CH 3 CCH lines 
from iMartfn et all (l2006bl) in the 2 mm band. For M 82, we 
complemen t our CS and CH 3 CCH detections with those of 
iBavet et all d2009l) and lMauersberger et all d!991l) . respectively. 
Finally, for IC 342 and Maffei 2 only few CS lines were taken 
from the literature (see Table [AT)- It has to be noted that some 
CS lines compiled from the literature were observed towards 
slightly different positions, so there might be a small error as- 
sociated to the different emission intensity in regions separated 
by a few arc seconds. In particular, the CS lines in NGC253 
differ a maximum of 6" in declination, while the CS lines in 
Maffei 2 differ a maximum of 7". In Sect. I3.ll we estimate the 
error introduced in our results by the difference in the observed 
positions. 

We found that for NGC 253 some of the compiled lines have 
a velocity resolution high enough for separating the two velocity 
components, roughly at 180 and 280 kms" 1 , which arise from 
the molecular l obes separated by 10", and located on both side s 
of the nucleus (Harriso n et al.ll 19991 [Mauersber ger et al.l l2003). 
For other observations, only one component was reported. In or- 
der to keep all data homogeneous, in our detections we have also 
fit only one Gaussian profile per line. In any case, from the two 
velocity comp onents analysis carried out bv lBavet et alj (2009) 
for CS and by IMartfn et all d2006bl) for HC 3 N, it can be seen 
that the differences between both components are of little signif- 
icance. 

In a similar way, Maffei 2 also shows two velocity compo- 
nents that could not be separated in the lines taken from the liter- 
ature. Therefore, in our detections we also fit only one Gaussian 
to both profiles. We have checked that a single Gaussian fit, 
both in NGC 253 and Maffei 2, leads to the same results as two 
Gaussians fit in terms of structure and properties of the gas. 
Thus, using just one velocity component does not affect our con- 
clusions. 



3. Results: Multi-transition analyses 

In this Sect, we present the detailed analysis of the multi- 
transition observations of CS, CH 3 CCH, and HC 3 N both un- 
der LTE and non-LTE approximations. During the analysis we 
derive the physical parameters for different molecular gas com- 
ponents within each source. It is important to note that these 
"gas components" are a simplification of the real structure of the 
molecular clouds. Though the real molecular clouds present gra- 
dients in temperature, density, and molecular abundances, this 
approach helps sketching the differences between the structure 
of molecular clouds in each galaxy. 



3.1. LTE analysis: Determination of rotational temperatures 
and column densities. 

We analyzed the detected transitions under the LTE approxi- 
mation, assuming optically thin emission. Estimations of the 



column density (AO and rotational temper ature (r rot ) were ob- 
tained using the Boltzmann diagrams (see Goldsmith & Langer 
1999 for a detailed explanation of the method and equa- 
tions). The spectroscopic parameters were taken from the 
CDMS dMuller et alJl200lllMuller et alj|2005l) and JPL catalogs 
dPickettet al.lll998l) . 

In the case where all lines are emitted under LTE condi- 
tions with an uniform excitation temperature, a single straight 
line can be fitted to all population levels. The derived 7^ will be 
a lower limit to the kinetic temperature (T^\ n ) if the lines are sub- 
thermally excited, e.g. if the H 2 densities are not sufficiently high 
to counterbalance spontaneous decay of the excited levels, which 
is the assumption used here. However, molecular emission arises 
from clouds with gradients in physical conditions which appears 
as different rotational temperatures components. In those cases, 
it is not possible to fit a single T lot , but several values of T m are 
needed to fit all the data. Although a continuous temperature dis- 
tribution seems more realistic, we here fit a discrete number of 
temperature components. This is sufficient to print out the tem- 
perature changes to be expected. In order to sample a wide range 
of physical conditions in the clouds, it is necessary to observe 
several transitions well separated in energies. This implies the 
observation of low and high excitation molecular lines. In this 
paper, we present the first detection of several high J HC 3 N and 
CH 3 CCH transitions, which allow us to sample for the first time 
the high excitation gas components previously undetected. 

The Boltzmann diagrams for all galaxies and molecules are 
shown in Fig. The column densities and rotational tempera- 
tures fitted through the LTE analysis are shown in Table [3] The 
column densities are source averaged with the source size of 
Table Q] Note that each N and T Iot are obtained from one lin- 
ear fit, so they refer to two or more population levels. We used 
the MASSA0 software package for this analysis. Relative abun- 
dances derived for the different gas components traced by every 
molecule are also shown in Table [3] 

As already pointed out in Sect. [2] the CS lines in NGC 253 
and Maffei 2 were observed towards slightly different positions. 
We have estimated the error introduced in the column densi- 
ties and rotational temperatures due to the differences in the ob- 
served positions. For NGC 253, CS lines differ a maximum of 
6" in declination. In the worst case, this represen ts 60% of the 
IRAM beam size. [Mauersber ger & Henkell (1 19891) observed the 
CS(2 - 1) line towards several positions of the NGC 253 cen- 
ter, separated by +10" both in declination and right ascension. 
Our (A, A)=(0",0") position corresponds to their (0", +3") po- 
sition in their map. Then, we can calculate the intensity change 
due to an offset of 10" in declination. Taking their (0",0") and 
(0", +10") positions, we obtain a decrease of 55% in the inte- 
grated intensity and 34% in the line intensity. Therefore, if the 
difference in declination is only of 6", one could expect that 
the integrated intensity, used for plotting the Boltzmann dia- 
grams, had to be decreased by less than 33%. If we re-plot the 
Boltzmann diagram using +33% of the integrated intensity of 
our CS(2 - 1), the variation of the column density and rotational 
temperature is less than 41% and 25% respectively. These vari- 
ations in Ncs and T rot also reflect the normal deviation of the 
resul ts between our work and that of Mauersberger & Henkell 
(119891 due to different instrumentation used and calibration er- 
rors. 

Unfortunately, we cannot make such an estimate for the 
CS(5 - 4) line, since it has not been observed towards several 
positions within the center of NGC 253. However, for a sepa- 



3 MASSA http://damir.iem.csic.eS/mediawiki-l.12.0/index.php/MASSA_Usei 





Fig. 5. Boltzmann plots obtained with MASSA for all the galaxies and molecules. The column densities and rotational temperatures 
of each component were obtained from the linear fits. The errors correspond to those of the integrated intensities of the Gaussian 
profiles. The errors in the CS diagrams of NGC 253 and Maffei 2 also take into account the uncertainty due to the different observed 
positions used for each line, as detailed in Sect. 13.11 



ration of 6" in declination, about half of the beam size at this 
frequency, we might expect a variation of less than a factor 1.3 
in column densities and 1.2 in rotational temperatures. 

In the case of Maffei 2, there are neither CS lines observed in 
several positions. The separation between the CS lines we have 
used is 7" and 4" in declination for the CS J=2 - 1 and 3 - 
2 respectively (1/4 of the beam in both cases, which means a 
variation of 15% of the integrated intensity). Thus, the expected 
difference in column densities and rotational temperatures are of 
less than a factor 1.2. 

The main results from the LTE analysis are: 

- CS and HC3N show similar excitation structure in NGC 253 
and IC 342, where three temperature components can be fit- 
ted. This fact suggests that these two molecules are well 
mixed and tracing the same gas components. This possibility 
is further discussed and supported by the non-LTE analysis 
carried out in Sect. 13.21 

- The excitation of the molecular gas in Maffei 2 and M 82 as 
traced by HC3N can be fitted by a single rotational tempera- 
ture. Only CS shows two different T Iot . The temperature gra- 
dient in these galaxies is significantly less pronounced than 
in NGC 253 and IC 342. Nevertheless, since we cannot rule 
out higher excitation gas components, we computed upper 
limit values of r, ot and column densities from the undetected 
highest J-transitions of HC3N and CH 3 CCH in Maffei 2 (see 
Table|AB. 



- CH3CCH shows a high excitation temperature relative to 
the other species for the galaxies in our sample (r rot = 13 
and 44 K for Maffei 2 and NGC 253 respectively, and 28 K 
for M 82). This is not conclusive for IC 342, since the er- 
rors are still too large. These high rotational temperatures 
are expected from t he relatively low d ipole moment of 
CH 3 CCH of 0.78 D dBurrelletal.| [T980h. compared to 3.7 
D for HC 3 N dDe Leon et al.lll985l) . which makes T rot from 
CH3CCH a better approximation to kinetic temperature than 
higher dipole molecules. Unlike other complex molecules 
such as CH3OH, methyl acetylene shows high column den- 
sities and relative abundances in M 82 (N = 8.5 x 10 14 cm" 2 , 
N/Nh^ — 1.1 x 10~ 8 ). Since this molecule mig ht be easily 
photo-dissociated by intense UV fields dFuente et al.i r2005). 
lower abundances would have been expected in a PDR nu- 
cleus like the center of M 82. 

- The high J lines of HC3N trace the warmest gas, reach- 
ing rather high rotational temperatures of 73. 3+14. OK in 
NGC253 and 70.6+21.2K in IC342. No other molecules 
have been found to have such high rotational temperatures 
in these galaxies so far, with the exception of those derived 
from ammonia of 120 K and 443 K for NGC 253 and IC 342, 
respectively. These results turn HC3N into one of the best 
suited species for tracing the warmest and densest molecular 
gas in galaxies. 
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Table 3. Parameters derived from the Boltzmann Diagrams. 



Galaxy 


Molecule 


N 
(cm- 2 ) 


* TOl 

(K) 


N/N H2 




NGC 253 




CS T 


(4.4+1.8) x 10 14 
(1.5+0.6) x 10 14 
(0.7 + 0.3) x 10 14 


4.4 ± 1.1 
12.0 ± 3 
33.1 ± 8.3 


1.1 x 10 


-a 




HC 3 N 


(8.1 ±2.9) x 10" 
(3.6+ 1.4) x 10 13 
(2.2 ±0.7) x 10 13 


11.6 ± 1.8 
26.1 ±4.5 
73.3 ± 14.0 


2.2 x 10 


-v 




CH 3 CCH 


(3.2 ± 1.0) x 10 14 


44.4 ± 7.7 


5.2 x 10- 


-v 


M82 




CS 


(2.0 ± 0.1) x 10 14 
(3.6 ± 0.5) x 10 13 


5.8 ±0.1 
15.1 ±0.9 


3.0 x 10 


->) 




HC 3 N 


(2.5 ± 1.1) x 10 l j 


24.7 ± 3.9 


3.2 x 10- 


11) 




CH 3 CCH 


(8.5 ±0.9) x 10 14 


28.1 ± 1.2 


1.1 x 10- 


-a 


IC 342 




CS 


(1.4 ±0.1) x 10 14 
(6.0 ± 0.1) x 10 13 
(2.3 ±0.1) x 10 13 


4.7 ± 0.1 
10.6 ± 0.2 
20.1 ±4.0 


3.8 x 10 


-9 




HC 3 N 


(2.7 ± 1.1) x 10 I J 
(1.6 ± 0.7) x 10 13 
(7.2 ± 3.2) x 10 12 


13.1 ±2.3 
20.6 ± 3.3 
70.6 ±21.2 


8.6 x 10~ 


1U 




CH 3 CCH * 


4.5 x 10 14 


70.0 


7.8 x 10- 


-9 


Maffei 2 




CS T 


(1.2 ± 1.0) x 10 14 
(6.9 ± 5.7) x 10 13 


5.4 ± 4.5 
9.7±8.1 


4.3 x 10 


-9 




HC 3 N 


(4.3 ± 0.6) x 10' j 
< 4.4 x 10 12 


11.6 ±0.8 
<53.1 


9.8 x 10- 


11) 




CH 3 CCH 


(9.0 ± 3.0) x 10 l j 
< 1.1 x 10 14 


12.8 ± 1.9 
<46.9 


2.0 x 10 





Column densities (N) and rotational temperatures (r iot ) obtained from the linear fits shown in Fig.[5]for each galaxy and molecule. 
The last column shows the relative abundances using the total H 2 column densities of TablefT] Before comparing, we have added up all the col- 
umn densities we have obtained for a molecule. 

The large errors in both N and r rot reflect the uncertainty in the worse case scenario due to the different observed positions for the CS lines 
(see Sect. ED- 

* The obtained values of r iot and A'cttjCCH are not reliable for the observed lines due to their large uncertainties (see the corresponding Boltzmann 
Diagram in Fig.[5}, so we have adopted T mt =70 K and calculated A'ch,cch through this value. 



3.2. Non-LTE analysis. Estimation of the H 2 densities. 

In order to quantify the volume densities, n(H2), and sup- 
port/explain the multiple rotational temperatures, we have used 
a Large Velocity Gradient (LVG) statistical equilibrium model. 
Excitation effects due to line radiation trapping were taken into 
account, using the LVG approximation consisting of molecular 
clouds with homogeneous spherical distribution. In this model, 
we had two free parameters: the volume density, n{H2>, and 
the column density per line width, N/Av. We assumed a back- 
ground temperature to 2.7 K, neglecting the contribution of any 
additional radiation field. Unfortunately, there are not collisional 
cross sections available for CH3CCH so we could not run our 
LVG model in that ca se. Therefore, for th i s mol ecule we have 
used a LVG model of Mauersberger et al.l (119911) . which is ex- 
plained in Sect. 13.2.21 For CS and HC3N, the collisional rates of 
our LVG model cover the first 12 energy levels for the first one 
and the first 23 energy levels for the second one. 

3.2.1. The kinetic temperature 

LVG results suffer to some extent from a degeneracy between 
the derived n(H{) and the kinetic temperature. Some molecules, 
like ammonia (NH3) or formaldehyde (H2CO), are outstand- 



ing species usually used as thermometers, since the relative 
populations of their metastable levels are essentially only sen- 
sitive to the 7\i n . Considering that the emission from ammo- 
nia and the molecules in our study arise from similar regions, 
the LVG degeneracy can be broken, since we can use the ad- 
ditional in formation on the T\m deriv ed from ammonia obser- 
vations by Mauersberge r et al.l (I2003I) . The assumption of am- 
monia arising for a similar volume that the other molecules is 
justified by the fac t that the large ammonia abundances found 
in these galaxies (Mauersberger et al. 2003) makes it very un- 
likely that a large fraction of the dense gas observed in CS and 
HC3N does not contain NH3. The exception is M82, where the 
ammonia abundance is rather lower than in the clouds of the 
Milky Way. In this galaxy, H2CO seem s to arise from a simi- 
lar volume than CO (Mti hle et al.l l2007). which indicates that its 
molecular emission is more extended than that of NH3, and also 
traces a warmer and more diffuse component in M 82 than NH3, 
CS, HC3N or CH3CCH. In fact, all these molecules have high 
dipole moments, need cri tical densities ~ 10 4 cm 3 to be excited 
(Mauersberger et al. 2003) much larger than those derived from 
H2CO dMiihle et al.ll2.007h . Thus, we find ammonia the most re- 
liable tracer of the kinetic temperature, since its characteristics 
(dipole moment, critical density and emission size) are in an bet- 
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Fig. 6. Results of the LVG model computations. The Y-axis is the log(nn 2 ) and the X-axis, the log(N), where N is the column 
density of the molecule. Each line represents the brightness temperature of a transition. When two or more lines cross, there is a 
density component. Elliptical shapes indicate where is more probable to find one of these density components within an uncertainty 

ofr b +io%K. 
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Table 4. Parameters obtained using our LVG code for CS and HC 3 N, and fhe lMauersberger et all (Il99ll) LVG model for CH 3 CCH. 



Galaxy 


Tkin 


Molecule 




n(H 2 ) 
(cm -3 ) 


N 
(cm ) 


^LTE / N non-LTE 


NGC253 


120 K 
















p c 


(1.7 


1 a £.\ \/ 1 a4 
± U.O) X IU 


(/.J ± 1. /) X iu 


A A 

U.O 








(3.1 


_i_ 1 A \ \/ 1 r\5 

± 1.4) X 1U 


/1 t _l n t\ v i a14 
{L.J + U.J) X IU 


1 1 

1 . 1 








(1.4 


1 A H \ * 1 A6 

± U. /) X IU 


10 a 1 1 a\ w 1 nil 

(o.U ± 2.0) X 11) 


0.9 






TIP TvT 


(2.0 


1 a Q\ w 1 A4 

± U.o) X 1U 


{ 1 .L ± 1 . 1) X IU 


1 1 
1 . 1 








(1.2 


± U.4) X IU 


(2.0 ± U.J j X 1U 


1 A 
1 .4 








O Q 

(2.y 


, 1 q\ w 1 a5 
± 1. j) X 1U 


1 . / ± U.J ) x iu 








PTI fYM I 

L.rl3L.Crl 




— IU 

1 n6 
~ 1U 






JVL o2 


3U K 


















(1.1 


± U.4) X 1U 


(J.Z ± 1.4) X IU 


U.O 








(2.3 


± 1.6) x 10 


(3.6 ± 1.3) x 10 


1.0 






HC 3 N 


(1.7 


± 1.3) x 10 


(1.9 + 0.6) x 10 


1.4 






PIT ppiT 




~ tu 






IC 342 


53 K 
















CS 


(3.3 


± 1. 4) X 10 


/l c 1 a n\ w 1 aJ4 

(2,5 ± U.9) X IU 


0.6 








(3.4 


± 1.6) X 10 


1 1 o\ \s 1 a13 
(0.2 ± 1.2) X lu 


1 A 
1 .U 








(1.2 


± 0.4) x 10 


/1 A . AC\«, f aH 

(1.9 ± 0.5) x 10 


1.2 






TIP T\T 


(0. / 


± Z.U) X 1U 


(OA xflAW I A1J 
(Z.4 ± U.O) X 1U 


1 1 
1 . 1 








(l.J 


i a c\ w i a5 
± U.J j X 1U 


(l.j ± U.j) X 1U 


1 . 1 








(2.S 


1 an v i a5 
+ U.O) X IU 


(y.i ± i.j) x iu 


U.cS 






PIT PPI I 

L.rl3L,L-rl 




~ IU 
~ 1U 






JVlalTCl 2 


A C If 
4o IV 
















PC 

Lb 


/ 1 Q 

(1.3 


_i_ A A \ \/ 1 A4 

± U.4) X 1U 


(S. O _i_ 1 ^\ w 1 Al4 

p. 2 ± l.jj X 1U 


A O 








(p. 5 


, 1 i \ v ia5 
± 1 . 1) X 1U 


(^4.0 ± U.V ) X 1U 


l.J 






HC 3 N 


(3.4 


± 2.9) x 10 3 


(1.5 ± 0.4) X 10 id 


2.9 






CH 3 CCH 




~ 10 4 






Maffei 2 


132 K 
















CS 


(4.5 


± 1.4) x 10 4 


(8.7 ± 2.8) x 10 14 


0.1 








(1.8 


± 0.6) x 10 5 


(4.6 + 0.9) x 10 13 


1.5 






HC3N 


(1.3 


± 1.1) x 10 5 


(1.5 ± 0.5) X 10 ]i 


2.9 






CH3CCH 











The last column reflects the ratio between the column densities obtained by both LTE and non-LTE approaches. 



ter agreement with CS, HC3N or CH3CCH than other tracers like 
CO and H 2 CO. 

Thus, we have used the kinetic temperatures derived by 
Mauersbe rger et al.l d2003l) as an input parameter in our LVG 
model. They obtained a r rot in NGC253 of 142 K and 100 K, 
for the two velocity components. Since we do not make any dis- 
tinction between these two velocity components, we have used 
the average value of rki n =120K. In the case of M82, the kinetic 
temperature traced by ammonia was 29 K. On the other hand, 
Maffei 2 and IC 342 show two different rotational temperatures, 
one for the lower (J, K) transitions of ammonia, and another one 
for the higher ( J, K) transitions. In the case of IC 342, the T rot 
are of 53 K and 443 K respectively. The first one is derived from 
the transitions (J, K)-(l, 1) to (3,3), having energies in their 
Boltzmann diagram ranging from to 150 K, while the other 
component of 443 K is obtained from the fit to the highest tran- 
sitions, (J, K) = (5, 5) to (9, 9), with energies from 300 K to more 
than 800 K . In our IC342 population diagrams, the energy levels 
reach, as maximum, 130K. Thus, one could guess that the gas 
component at 443 K has little contribution to the relatively low 
energy transitions we deal with in this paper. So, we have used 
2\in= 53 K in our LVG analysis shown in Table[4] 

For Maffei 2, the two 7\j n of 48 K and 132 K were taken into 
account when running the LVG model. Table [4] shows the LVG 



results for both temperatures. Gas densities obtained from CS 
for this source are more sensitive to the changes of the 7\i n than 
those obtained from HC3N, with variations of up to one order 
of magnitude in the lower density component. Moreover, the CS 
and HC3N transitions seem to trace gas with low excitation tem- 
perature, as shown by the LTE analysis in Sec 13.11 (from 5 to 
10 K for CS and 12 K for HC3N). Following the same arguments 
used for IC 342 when comparin g our upper energy values o f the 
Boltzmann plots with those of Mauersberg er et all (|2()03). and 
taking into account that the rotational temperatures are consid- 
ered as a lower limit of the kinetic temperature, it seems more 
appropriate to use lion =48 K instead of Tki n = 1 32 K for Maffei 2. 



3.2.2. LVG model results 

Fig. [6] shows the predictions from our LVG modeling for all 
galaxies in a log[n(H2)] - log[N] diagram for CS or HC3N and 
taking the kinetic temperatures given above. Lines represent the 
model results matching the observed brightness temperatures for 
each transition. The points where two or more lines intersect 
correspond to the H2 density and the column density which fit 
those lines. The dashed ellipses show the area where two or 
more lines intersect to within an error of ±10%. Such conser- 
vative error was chosen taking into account the T^ n errors ob- 
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tained by Mauersberger et al.l (120031) . which are always lower 
than 10% (with the only exception of the 30% uncertainty in the 
Trot =48 K for Maffei2. Nevertheless, such error does not sig- 
nificantly affect our derived parameters). As shown in Fig. 6, we 
find that several components with different de nsities are needed 
to fit all observed lin es (as already shown bv lBavet et al.ll2008l 
and lBavet et al.ll2009l) . in agreement with the multiple rotational 
temperatures obtained from the Boltzmann diagrams. Table [4] 
shows the derived nfHz) and N for the different density compo- 
nents in all the galaxies. The errors in H2 densities and column 
densities are derived from the lcr errors in the Gaussian fits. The 
ratios between the column densities of CS and HC3N derived 
from the LTE analysis and the LVG model for the different den- 
sity components are also shown in the last column of Table [4] 

Unfortunately, our LVG model does not include the colli- 
sional cross sections of CH3CCH, and we do not have at our 
disposal any other similar code prepared to use it. So far, the 
best ap proximation to mod e l the excitation of CH3CCH was 
done bv lMauersberger et alj dl991f) . They run a LVG model for 
methyl acetylene using the CH3CN collisional rates and rota- 
tional constants, since these parameters may be similar for both 
sp ecies. Therefore, fo r meth yl acetylene we used the results 
of Mauersberg er et ail (Il99lh as a rough estimation of the H2 
volume densities within an order of magnitude uncertainty, i.e. 
~ 10 4 ,~ 10 5 or ~ 10 6 c m~ 3 . On the other hand, al though the 
temperatures used in the Mauersberger et afj dl99ll) model for 
CH3CCH do not exactly match those assumed from NH3, these 
differences in the 7\i„ do not play a relevant role and the vol- 
ume densities hardly vary when using either the ammonia or the 
methyl acetylene kinetic temperatures. 

The main results from both LVG analyses can be summa- 
rized as follows: 



- LVG analyses have been used to quantify the density and the 
molecular column densities of the gas components inferred 
from the Boltzmann plots. The column densities obtained by 
both approaches are consistent within a factor of 2 in almost 
all cases (see TableHJ). 

- CS traces a very wide range of densities. Similar to the 
LTE results our data requires three density components in 
NGC 253 and IC 342, with a density contrast of up to two or- 
ders of magnitude between the more diffuse gas (~ 10 4 cm~ 3 ) 
and the densest material (~ 10 6 cm~ 3 ). On the other hand, 
only two density components are needed to fit the data in 
Maffei 2 and M 82, with a smoother density gradient of about 
one order of magnitude. 

- HC3N traces a narrower range of densities than CS in M 82 
and Maffei 2. In these sources, only one density component 
is found, with a density of a few 10 6 cm 4 in M82, and 
a several 10 5 crrT 3 in Maffei 2. In fact, cyanoacetylene is 
not found in the more diffuse gas (~10 4 cm 4 ) observed in 
CS in these galaxies. For the other two galaxies, IC 342 and 
NGC 253, both species trace the same density components, 
but HC3N systematically provides lower densities than those 
ofCS. 

- Using the grid of models of Mauersber ger et all (120031) . we 
have found that CH3CCH seems to be arising from an inter- 
mediate component of ~ 10 5 cm 4 in M 82, while for IC 342 
and NGC 253 could be tracing gas of densities ~ 10 5 cm 4 
or/and ~ 10 6 cm 4 . In the case of Maffei 2, methyl acetylene 
arises from less dense molecular gas of densities ~ 10 4 cm -3 . 



4. Discussion: Molecular cloud structure in the 
nuclei of galaxies 

In order to understand the physical structure and properties of 
the ISM of the galaxies in this paper, we will use the simple 
molecular cloud structures sketched in Fig. [7] It illustrates the 
expected evolution of molecular clouds in starburst galaxies. We 
consider three different regions: The diffuse outer region (halo) 
with densities of ~ 10 4 cm 4 referred to as Region I; the dense 
gas with densities of ~ 10 5 cm 4 referred to as Region II; and the 
inner and densest gas we can trace (core), the Region III, with 
densities of ~ 10 6 cm 4 . Fig. illustrates how the CS, HC 3 N 
and CH3CCH seem to be distributed in the molecular clouds of 
the central hundred parsecs of the galaxies in our sample. At 
first glance, the molecular clouds of NGC 253 and IC 342 have a 
very similar density distributions, while Maffei 2 and M 82 show 
a different structures. Next, we will discuss why these density 
distributions point to three different stages of the starburst evo- 
lution. 

We can make a rough estimation of the molecular cloud 
sizes, which allows us to know how many molecular clouds 
could be contained in our beam. Supposing that the gas is in 
a virial equilibrium and formed by homogeneous clouds of uni- 
form H2 density of n = 10 5 cm 4 , we have obtained an approx- 
imate molecular clouds radii of ~29, ~51, ~19 and ~34pc for 
M82, NGC 253, IC342 and Maffei 2 respectively. Taking into 
account the source size used for each galaxy (see Table [TJ, it is 
possible to estimate the number of molecular clouds contained 
in the beam of about ~7, ~7, ~9 and ~5 for M82, NGC 253, 
IC342 and Maffei 2. Although there are not big differences 
in the numbers, they seem to be consistent with the proposed 
overall scenario: M82, the galaxy most dominated by the UV 
fields, is claimed to have smaller and more fragmented clouds 
(lFuenteetalJl2008h . Maffei2, on the other hand is the galaxy 
most dominated by shocks and its clouds seem to be rather big- 
ger, having less number of them in the beam (~5). However, 
these are just rough calculations, since we show in this paper the 
non negligible temperature and density gradients within molec- 
ular clouds, as depicted in Fig. [7] In fact, the observed gas could 
well be formed of several molecular clouds of different sizes 
each. Thus, this is something that could be better constrained us- 
ing theoretical models or interferometric data. Also, the assumed 
distances to the galaxies play an important role in the values ob- 
tained. Nevertheless, in order to see whether our values are con- 
sistent with those of the literature, we have compared the most 
straightforward case of IC342. The size of the molecular clouds 
we have obtain ed in this galaxy and the HC3N(10 - 9) interfer- 
ometric map of iMeier & Turner! (120051) are in good agreement, 
since they point out that the emission of HC3N in IC 342 is not 
expected to be more extended than ~50 pc. 

On the other hand, from the analysis presented in Sect. [3] 
the CS abundances remain almost constant independently of 
the dominant process, with the exception of NGC 253. This 
is likely due to the fact that CS is not significantly affected 
by any of the dominant p rocess (PDRs or shocks), as shown 
bv iReauena- Torres et al.l (1200 6) in the Galactic Center, and 
iMartfn et a l. (2009a) in a sample of galaxies. On the other hand, 
HC3N shows low abundance in M 82. This molecule seems to 
be easily dissociated in PDR s like in the center of this ga laxy, 
as was already proposed by Rodrigue z-Franco et al.l (119981) . At 
the same time, HC3N is absent in the non-detected Region I of 
M 82, likely due to the lack of shielding from the dissociating 
radiation. This is consistent with the idea that CS is more resis- 
tant to the UV than HC3N; and cyanoacetylene was not detected 
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Maffei2 N6C253 & IC342 M82 




Evolution of Starbursts 



Fig. 7. This simple cartoon shows how CS, HC3N and CH3CCH could be distributed in the giant molecular clouds within the central 
hundred parsecs of the four galaxies from the LVG model results. The size scales range from ~145 to ~560pc depending on the 
galaxy distance and beam size. We have found a Region I of densities about ~10 4 cirT 3 , an intermediate Region II of densities 
~10 5 cirT 3 , and a Region III of high densities around ~10 6 cirT 3 . Maffei2 shows two density components, but no core (Region III) is 
detected. In NGC 253 and IC 342, up to three gas components of different densities are distinguished (~10 4 , ~10 5 and ~10 6 cm -3 ). 
M 82 shows a region of intermediate densities and a core. The Region I is not detected by the molecules used here, although it 
should be surrounding the Region II. 



Table 5. N(HC3N)/N(CS) ratios and percentage contribution of each Region to the overall dense molecular gas 





Maffei2 


IC 342 


NGC 253 


M82 






(A'hcn/A'cs - 


%contribution) 




Region I (~10 4 cm^) 


... - 92% 


0.10-75% 


0.10-77% 




Region II (~10 5 cnT 3 ) 


0.33 - 8% 


0.39 - 19% 


0.33 - 14% 


...-90% 


Region III (~10 6 cm- 3 ) 




... - 6% 


... - 9% 


0.53 - 10% 



The fist value indicates the N(HC3N)/N(CS) ratio obtained from our LVG code. The second value represents the percentage contribution of each 
Region (I, II or III) to the overall dense molecular gas, traced by CS, within the molecular clouds in the center of the galaxies, obtained from 
our LVG model. 



in the Region I in Maffei 2 either. In this case, the most likely 
explanation is that the low density in this region does not excite 
the molecule, so its lack could be due to excitation effects rather 
than to chemistry. 



As shown in Table [5] the HC3N/CS ratio hardly varies when 
comparing regions with similar densities among the galaxies. 
This supports our sketch of the molecular clouds divided in re- 
gions, as they not only share similar densities but also molecular 
abundances. In Table [5] we also present the fraction of gas that 
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each of these components contribute to the overall molecular gas 
traced by CS (we take this molecule as a reference since it is the 
one which traces a wider range of densities). In the following 
we discuss the structure of the different observed sources and its 
possible relation to the evolution of their nuclear starburst. 

4.1. Molecular clouds in galactic centers hosting young 
starburst. The case ofMaffei2 

Maffei 2 is considered to be at a very early state of the starburst 
evolution (Ma rtrn et al.ll2009ah . with a large amount of reser- 
voir gas in its center, and only a moderate star formation rate 
(Table Q]). The gas kinetic t emperatures derived from ammo- 
nia ( 7\i n = 48 and 1 32 K from lMauersberger et al .1120031) indicate 
that shocks dominate th e ISM heating (as already indicated by 
e.g. llshiguro et al.ll9 89). which is consistent with the HNCO/CS 
ratio. There are several processes that can contribute to generate 
shocks in the early phases of the starbursts. One is the existence 
of a bar, which creates a potential leading to cloud-cloud col- 
lisions at the edges. This is, in fact, the case of Maffei 2. The 
molecular cloud collisions can be also generated by the usual 
compression which leads to the star formation; and there could 
also be a contribution from the stellar winds of the newly formed 
stars. 

The molecular clouds in Maffei 2 seem to have large ha- 
los (i.e., Region I, which contains > 90% of the detected CS) 
of molecular gas very rich of fragile molecules, like for exam- 
ple am monia or HNCO, but also CO, as seen by Israel & Baas 
(2003). The molecular complexes of Maffei 2 show low rota- 
tional temperatures of 5-15 K in the molecules studied, and 
normal fractional abundances of ~ 10~ 9 for CS, HC3N and 
CH3CCH (see Table O, which are similar to those of the "qui- 
escent" molecular clouds in the Galactic Center (GC), with no, 
or low, star formation activity. On the other hand, if we consider 
the upper limits to the undetected lines in Maffei 2 in the LTE 
analysis, the derived rotational temperatures are less than 53 K 
from HC 3 N and 47 K from CH3CCH (see TableEJ. 

Taking into account all the results from the LTE and non- 
LTE analyses, the three studied molecules seem to be arising 
from molecular clouds similar to the molecular clouds with low 
star formation activity found in the central region of the Milky 
Way. At this stage, the molecular material is likely to be fun- 
neled towards the center of the potential well. In fact, a bar 
of molecular gas drives the material to the inner central parts 
dlshiguro et all 19891 iHurt & Turner! 199ltlKuno et alJ2008l) . in- 
creasing its density and feeding the still relatively low star for- 
mation rate found in Maffei 2. 



4.2. Molecular clouds in galaxies hosting intermediate-age 
starbursts: NGC253 and IC342 

The two galaxies in our sample, NGC 253 and IC 342, which are 
considered to host a starburst at intermediate stage of evolution, 
where a large amount of gas has been already converted into 
starts dRieke et all 119881: ISchinnerer etail 120081: iMartm et al.1 
2009bj), show molecular clouds with similar physical and chem- 
ical structure (see Table |5J. The median age of the star clus- 
ters in th e center of NGC253, studied at o ptical wavelengths, is 
~6Myr dFernandez-Ontiveros et aL 2009), while fro m near in- 
frared (NIR) continuum emission, Boke r et al.l (119971) found that 
the center of IC 342 seems to be dominated by a young cluster of 
~10Myr. Although the optical and the NIR ranges trace stellar 
populations with different ages, these estimations seem to indi- 



cate that the star bursts in both sources have probably started at 
a similar time ago. 

Both NGC 253 and IC342 show a very high density con- 
trast with nearly two orders of magnitude between the Region 
I and the Region III. Furthermore, Regions III are rather small 
as compared with the Regions I, but are likely responsible for 
feeding the star formation at high rates in these galaxies. On the 
other hand, CO arises in this galaxies from gas w i th densities 
in the range 10 3 - 10 4 cm" 3 dlsrael & Baasl 120031: iBavet et al.l 
12004 I2006t lGustenetal]|2006l) . The average clouds in these 
two galaxies resemble the Sgr B2 star forming region in the 
GC. This cloud is centrally peaked sur rounded by a large molec- 
ular envelope (Ide Vicente et al.Nl 9961) and shows that very re- 
cent star formation is still going-on in the very high density 
core (Ide Vicente et al . 2000). As we have mentioned in Sect. [T] 
IC 342 shows evolved star formation in the center, while at l arger 
radii the stars seem to be younger dMeier & Tu rner 200J). We 
would like to emphasize that with single dish telescopes the 
overall scenario is a mixture of both cloud structures, where 
shock-dominated and PDR-dominated regions are not distin- 
guished. Thus, the resul t is a mixture that l eads to an interme- 
diate stage of evolution (Martin et al. 2009b), where still a large 
fraction of the gas is in a high density component and the cores 
are actively formin g stars. As suggested by the detection of PDR 
tracers in NGC 253 dMartm et al.l2009bl) and illustrated in Fig. El 
photodissociation starts playing a role in these clouds in their 
outermost layers. 



4.3. Evolved starburst: M82 

The stellar clusters in the center of M 82 are estimated to be the 
result of starbursts from a t least ~ 10- 15Myr ago, a s seen 
at optical wavelengths by Konstantopoulo s et al.l (|2009). This 
value indicates that we are seeing a more evolved starburst (or 
post-starburst) in M 82, older than the NGC 253 and IC 342 ones 
(we are not aware of a published stellar population age in the 
center of Maffei 2). On the other hand, it is well established 
that newly formed stars in M 82 create large ph otodissociation 
regions in its central region (Mar tin et al.ll2006bl) . which domi- 
nate over the other physical processes. Here, we follow the pro- 
posed overall scenario of M 82 as a giant PDR, although we 
are developing for the future a more sophisticated interpreta- 
tion where PDRs, XDRs and dense gas chemistry are mixed and 
will be better disentangle. This galaxy shows clouds with less 
density gradients than those found in the other galaxies, with 
a density contrast of only a factor of 6 between the two de- 
tected gas components (Regions II and III). This could be due 
to the fact that, in this evolved starburst galaxy, the H2 density 
in Region II, of about ~10 5 cm -3 , is surrounded by an enve- 
lope (not detected) where some molecules, especially the com- 
plex ones, are destroyed. The low density envelope surrounding 
the Region II could be dominated by mainly atomic composi- 
tion wit h small column de nsities of resistant molecules like CO. 
In fact, iMao et al.l d2000t) have detected CO in M 82 with den- 
sities of ~10 3 ' 7 and ~10 3 cirT 3 , suggesting that this molecule 
does not trace the same material than high density tracers like 
CS. Likely the mechanical energy injected by the evolved star- 
bursts had dispersed the outer parts of the clouds making larger 
and more diffuse envelopes. Thus, the number of gas compo- 
nents with different densities and temperatures can be under- 
stood within the starburst evolutionary framework; in the later 
stage the dense molecular tracers are not arising from the outer 
parts, since they are dissociated in such regions. At this stage, the 
outermost photo-dissociated envelope, of densities < 10 5 cm" 3 , 
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might represent a large fraction of the cloud, but with low abun- 
dance of complex molecules. At the same time, there is an im- 
portant contribution of the Region III (~10 6 cm 3 ), directly re- 
lated with the large amount of star formation that have already 
been taken place. The hi gh density compone nt is also shown by 
observations of CH 3 OH dMartm et ai]|2006bl) . 

Methyl acetylene (CH3CCH) shows a very different behav- 
ior in M 82 than the other molecules. First, its column density 
in M 82, of 8.5 x 10 14 cm" 2 , is a factor of ~4 larger than that of 
CS, and more than one order of magnitude than the HC3N one. 
Its high fractional abundance with respect to H2, 1.1 x 10~ 8 , is 
puzzling for a molecule which is supposed to be easily disso- 
ciated in PDRs dFuente e t al. 2005). The comparison with other 
molecular clouds dominated both by UV radiation and shocks, 
in the Galactic Center and in external galaxies, shows that M 82 
is the galaxy with the highest CH3CCH abundance observed up 
to now. SgrA* (-30,-30) and G+0. 18-0.04, used as templates of 
molecular clouds dominated by phot odissociating radiat ion in 
the Galactic Center of the Milky Way (Mart in et al.l2008al) . have 
relative abundances larger than a few 10~ 9 (Martin, 2006). On 
the other hand, the relative abundances we find in NGC 253 and 
Maffei2 are 5.2 x 10~ 9 and 2.0 x 10 9 respectively. One possi- 
ble explanation for the high CH3CCH abundance in M 82 is that 
this molecule can be created through gas phase ion-molecule or 
neutral-neutral reactions, while other molecules, like methanol 
or ethanol, are possib ly only formed on ice-layers of dust grains 
dBisschop et alj|2007h . When these complex molecules are dis- 
sociated by the strong UV fields of M 82, their abundances de- 
crease, while a similar decrease for CH3CCH is balanced by ion 
molecule reactions favoring its formation. On the other hand, a 
strong anti correlation bet ween Ncmccvi a nd 7\; n has been ob- 
served in galactic sources dLee et alTl9 96), and the gas in M82 
is indeed as cold as ~30 K. Methyl acetylene is not observed in 
the regions of densities < 10 5 cirT 3 since it would be dissociated, 
but in regions of intermediate densities (Region II), where it is 
well shielded from the dissociating radiation. Furthermore, we 
do not see it in the Region I, presumably because its abundance 
in such places is almost negl igible, although it could exist there 
dChrurchwell & Hollisll 19831) . 

5. Conclusions 

From our multi-line studies of CS, HC3N and CH3CCH towards 
the centers of the starburst galaxies, M 82, NGC 253, IC 342 and 
Maffei 2, and using both LTE and non-LTE approaches, we con- 
clude that: 

- CS and HC3N are tracing gas components with similar phys- 
ical properties. CS arises from a slightly denser gas (see 
Table|4]i, while HC3N arises from warmer gas (see TableO. 
Both molecules are excellent tools to derive the density 
structure of the average molecular clouds in galaxies. 

- CH3CCH arises from a "warm" (T mt from 13 K up to 44 K) 
gas components without large rotational temperature gradi- 
ents. This molecule shows very high abundances and column 
densities in M82 (1.1 x 10~ 8 cnT 2 and 8.5 x 10 14 cnT 2 re- 
spectively), and probably arises from a region of densities 
~10 5 cirT 3 which is well shielded from the strong UV radia- 
tion. 

- The gas traced by the observed molecules show differences 
between the galaxies in terms of density structure and com- 
position. This fact allowed us to link the evolution of the star- 
burst with the structure of the molecular gas clouds. First, 
Maffei 2, the galaxy with the youngest starburst, shows a 



large low density halo (Region I ~10 4 crrT 3 ) and a more 
denser Region II with densities of ~10 5 crrT 3 . This galaxy 
does not show signs of a dense gas core (Region III), prob- 
ably because its starburst is still starting. Then, in an inter- 
mediate stage of evolution, NGC253 and IC342 clearly show 
three different gas components. First, a halo of low densities 
(Region I~10 4 crrT 3 ), possibly surrounded by a small enve- 
lope where the molecules are dissociated. Then, an interme- 
diate zone (Region II) of densities ~ 1 5 cm" 3 ; and finally a 
core of very dense gas (Region III~10 6 cm 3 ) which points 
out to the on-going star formation. Finally, M 82, the more 
evolved starburst galaxy, shows a Region II of intermediate 
densities of ~10 5 cirT 3 , probably surroundedby a large enve- 
lope dominated by photodissociating radiation, and a Region 
III of denser gas, comparable in size to the cores of NGC 253 
and IC 342. 
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Appendix A: Table A.1. Gaussian fits parameters results. 
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HC 3 N(17- 16) 


3.4+0.5 


63.0* 


184.0 * 


50.0 






d 




(0,-9) 


HC 3 N(26-25) 


3.2+0.7 


63.0* 


184.0* 


47.0 






d 




(0,0) 


CH 3 CCH (8 - 7 ) 


3.4+0.1 












c 




(0,0) 


CH 3 CCH(9 - 8 ) 


4.3+0.2 












c 




(0,0) 


CH 3 CCH(10o-9 ) 


5.1+0.5 












c 




(0,0) 


CH 3 CCH(13 - 12o) f 


2.9+0.3 


155.0 


250.0 


17.3 


3.8 


90.0 


a 


M82 




(+13, +7.5) 


CS (2 - 1) 


9.4+0.2 


102.6 


300.1 


86.3 


4.6 




b 




(+13, +7.5) 


CS (3 - 2) 


11.5+0.2 


105.6 


293.8 


103.7 


2.6 


481 


a 




(+13, +7.5) 


CS (4 - 3) 


6.5+0.8 


87.7 


309.2 


70.1 


16.6 




b 




(+13, +7.5) 


CS (5 - 4) 


4.2+0.2 


71.6 


316.4 


55.5 


7.3 


248 


a 




(+13, +7.5) 


CS (7 - 6) f 


0.9+0.1 


53.0 


317.8 


16.7 






b 




(+13, +7.5) 


HC 3 N(15- 14) 


1.2+0.2 


97.3 


308.8 


11.5 


2.5 


59 


a 




(+13, +7.5) 


HC 3 N(16- 15) 


1.0+0.1 


100.0 * 


300.0 * 


9.6 


1.4 


292 


a 




(+13, +7.5) 


HC 3 N(17- 16) 


1.3+0.1 


96.9 


314.3 


13.0 


2.4 


158 


a 




(+13, +7.5) 


HC 3 N(18- 17) 


0.8+0.1 


87.5 


308.5 


8.6 


2.0 


217 


a 




(+13, +7.5) 


HC 3 N(24-23) 


0.4+0.3 


42.6 


333.8 


9.9 


4.6 


294 


a 




(+13, +7.5) 


HC 3 N (28 - 27) 


<0.2 


100.0 * 






1.1 


1400 


a 




(+13, +7.5) 


CH 3 CCH(6 -5 ) 


3.7+0.1 


123.8 


315.0 


28.1 


1.0 


168 


h 




(+13, +7.5) 


CH 3 CCH(8 -7 ) 


7.4+0.6 


122.9 


310.0* 


56.8 


2.5 


59 


a 




(+13, +7.5) 


CH 3 CCH(9 - 8 ) 


7.1+0.1 


109.6 


322.6 


63.9 


1.7 


158 


a 




(+13, +7.5) 


CH 3 CCH(10o-9 ) 


8.7+0.2 


97.5 


325.0 


84.1 


3.0 


281 


a 




(+13, +7.5) 


CH 3 CCH(15 - 14 ) 


7.4+0.7 


85.1 


339.1 


81.2 


8.0 


56 


a 




(+13, +7.5) 


CH 3 CCH(16 - 15 ) 


2.4+0.1 


70.1 


337.5 


32.1 


3.1 


799 


a 


IC342 




(0,0) 


CS (2-1) 


5.0+0.1 


53.9 


31.7 


87.4 


2.5 


28 


a 




(0,0) 


CS (3 - 2) 


4.9+0.1 


51.7 


31.2 


88.5 


3.0 


28 


a 




(0,0) 


CS (5 - 4) 


3.4+0.1 


63.4 


38.2 


37.6 


2.9 


57 


a 




(0,0) 


CS (7 - 6) t 


1.0+0.4 


75.6 


29.3 


12.3 


8.5 


50 


b 




(0,0) 


HC 3 N (9 - 8) 


0.7+0.1 


54.4 


34.4 


12.2 


1.2 


126 


a 




(0,0) 


HC 3 N(10-9) 


0.7+0.1 


40.9 


35.5 


15.3 


2.8 


84 


a 




(0,0) 


HC 3 N(11 - 10) 


0.7+0.1 


58.7 


34.0 


12.0 


1.0 


218 


a 




(0,0) 


HC 3 N(12- 11) 


0.7+0.1 


52.7 


39.3 


12.4 


1.3 


77 


a 




(0,0) 


HC 3 N(15- 14) 


0.5+0.1 


50.9 


39.2 


9.8 


2.0 


56 


a 




(0,0) 


HC 3 N(16- 15) 


0.4+0.1 


50.0* 


31.0* 


6.7 


1.8 


190 


a 




(0,0) 


HC 3 N(17- 16) 


0.7+0.1 


101.3 


44.1 


6.6 


1.5 


416 


a 




(0,0) 


HC 3 N (24 - 23) 


0.7+0.2 


76.7 


12.7 


8.4 


2.8 


211 


a 




(0,0) 


CH 3 CCH (5 - 4 ) 


0.2+0.1 


36.8 


33.0 


3.9 


5.8 


278 


a 




(0,0) 


CH 3 CCH(6o-5 ) 


0.3+0.1 


60.0* 


38.6 


4.5 


6.7 


276 


a 




(0,0) 


CH 3 CCH(8o-7 ) 


0.7+0.1 


50.0* 


25.5 


13.7 


1.5 


56 


a 


Maffei2 




(0,-7) 


CS (2 - 1) m 


0.9+0.1 


132.9 


-19.0 








f 




(0,-4) 


CS (3 - 2) 


6.1+0.5 


175.9 


-36.0 


33 






g 




(0,0) 


CS (5 - 4) 


3.0+0.3 


139.2 


-50.0 * 


20.3 


3.1 


99 


a 




(0,0) 


HC 3 N(9-8) 


1.0+0.1 


130.0* 


-33.3 


7.6 


78.1 


112 


a 




(0,0) 


HC 3 N(12- 11) 


1.8+0.4 


112.2 


-44.6 


15.3 


3.9 


21 


a 




(0,0) 


HC 3 N(15- 14) 


0.4+0.1 


130.0* 


-50.0 * 


3.1 


1.0 


112 


a 
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(0,0) 


HC 3 N (23 - 22) 


<0.3 


130.0* 






2.2 


183 


a 


(0,0) 


HC 3 N (28 - 27) 


<0.2 


130.0* 


_ 


_ 


8.5 


140 


a 


(0,0) 


CH 3 CCH (5 - 4 ) 


0.3±0.1 


142.5 


-50.0 * 


2.1 


1.0 


224 


a 


(0,0) 


CH 3 CCH(6 -5 ) 


0.8±0.1 


146.3 


-50.0 * 


4.8 


1.0 


224 


a 


(0,0) 


CH 3 CCH(8 -7 ) 


0.7±0.1 


130.0* 


-60.4 


5.3 


1.1 


112 


a 


(0,0) 


CH 3 CCH(13 - 12 ) 


<0.3 


130.0* 






3.0 


224 


a 


(0,0) 


CH 3 CCH(14 - 13 ) 


<0.2 


130.0* 






2.0 


253 


a 



Table A.l. All lines were observed with IRAM 30-m telescope except: ": line 
observed with SEST; f. line observed with JCMT; ®, line observed with 12- 
m NRAO. A * following the FWHM or (/and) V LSR means that the parameter 
was fixed during the Gaussian fit. Lines with values of only rms and integra- 
tion time are not detect ions, but upper limit s. Ref erences: (a): this work; (b) : 
iBavet et al.l <2009h: (c): iMartfn et al] d2006bl): (d): IMauersberger et al.l JT 990): 
(e): [M auersber ger et al.l d!991h : (f): ISage et al.l dl990l) : (g): IMauersberger et al.l 
<1989h : (h):Aladro et aL 2009b, in prep. 



